Realistic Derivation of Heisenberg Dynamics 

Paul J. Werbos* 
National Science Foundation^, Room 675 
Arlington, VA 22203 

ABSTRACT 

Einstein conjectured long ago that much of quantum mechanics might be derived as a 
statistical formalism describing the dynamics of classical systems. Bell's Theorem 
experiments have ruled out complete equivalence between quantum field theory (QFT) 
and classical field theory (CFT), but an equivalence between dynamics is not only 
possible but provable in simple bosonic systems. Future extensions of these results might 
possibly be useful in developing provably finite variations of the standard model of 
physics. 
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I. INTRODUCTION AND SUMMARY 



The goal of this letter is to present a new formalism for analyzing the statistical dynamics 
of "classical fields" - of ODE or PDE derived from a classical Lagrangian density of the 
usual Hamiltonian form. To keep the notation simple, I will focus on the ODE case with: 

n 

L = ±Y,(Vj n j- 7i j<Pj)- H (<P>x) (U 



where the state of the system at time t is defined by the two classical vectors $ , 7C g R n . 
The generalization to PDE is discussed in Appendix B. 

The following section will show how a probability distribution or statistical 
ensemble of possible states, Pr((|>, 7t) may be "encoded" or "mapped" into a "classical 
density matrix" p, which is dual to the usual field operators of the Heisenberg 
formulation of QFT 1 . 1 will then define simple field operators <£j and rij which have the 
properties that Tr(<£jP)=<(Pj>, Tr(]ljP)=<7lj> and both obey the usual 
Heisenberg/Liouville dynamics as described by Weinberg 1 . Note that I use angle brackets 
to denote the expectation value of classical stochastic variables. 

The concluding section will discuss the possible significance and future extension 
of these results. 

II. CONCEPTS AND DEFINITIONS 

Efforts to understand the statistics of classical fields have a long history which cannot be 
reviewed comprehensively in a brief letter; however, some historical concepts are 
essential to understanding of the new work. 

For centuries, statisticians have studied how probability distributions Pr((|>) can be 
represented and studied in terms of statistical moments, correlations or cumulants. For 
example, the set of real numbers 



includes the statistical mean and (implicitly) the covariance, skewness and kurtosis, etc. 
The set of all such u values for a given probability distribution may be seen, 




(2) 
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mathematically, as a kind of a real wave function, as a real vector u in the usual Fock- 
Hilbert space of QFT! 

Many classical physicists have tried to "close turbulence" by deriving linear or 
nonlinear dynamical equations for the statistical moments. However, when the vector u is 
used directly - without scaling and without exploiting its special properties - the 
dynamical equations tend to display a problem of infinite regress, explicitly or implicitly. 
Many of those efforts are analogous to older textbooks which try to solve simultaneous 
linear equations by writing each equation separately and manipulating them one by one; 
by shifting over to Fock space concepts from QFT, we can achieve a simplicity and 
power analogous to the power which matrix algebra provides in working with systems of 
linear equations. 

In a previous paper, we 2 suggested that the statistical moments are slightly easier 
to analyze if we first scale them: 



_ '■ .'2 ■■■■■'„ ( 2) 
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which is equivalent to the more compact and modern definition: 



v((p) = exp 



f n \ 

v '■=' J 



I0> , (4) 



where a/ is the usual creation operator of QFT. This lets us represent a probability 
distribution Pr($) by: 

v=<v(9)>=|v(9)Pr(£)J> (5) 
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Equation 4 is very similar to the mappings used by Segal and by Rajeev 3 in the analysis 



of Bose-Fermi equivalences; however, it is not a special case, because the operator in the 
exponential is not antiHermitian, 

It is easy to calculate Ivl 2 and verify that it is bounded but not normalized to 1 for 
most states of physical interest. Furthermore, v only represents a distribution Pr((|>), not a 
joint distribution Pr((|>,7t). Therefore, to represent the statistical moments of the system 
given in equation 1, 1 propose the following normalized (complex) vector: 



where the superscript H denotes the usual Hermitian conjugate. Vectors w which 
represent a specific state ($,7t) as in equation 6 will be called "pure states." Matrices p 
which represent an ensemble of states as in equation 7 will be called "physically 
realizable." 

For this ODE system, the two basic field operators take on a simple form: 




(6) 



The classical density matrix is then defined by 




(7) 




(8) 



n, 



i 



(a j - a)) 



(9) 
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III. PROPERTIES AND DYNAMICS 
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It is easy to verify that <£j and 11 possess the usual properties for field operators, apart 
from minor scaling factors, such as: 

^.,nj=^(})(i)(-2) = i^ (io) 



Likewise: 

b 7 ,(nj m J=|^m(nj 



m— 1 



(11) 



For any analytic function of $ and 7t, 



we may define two different "quantized" versions, a canonical version 



/ c (^,0)= £ c; : ; ci>:ci>v...*;:n;...n; 



and a normal version 



/»(*,n)= :/ c (^,n) 



(12) 



(13) 



(14) 



where the colons are the usual notation for normal product. It is straightforward but 
tedious to prove the following extensions of equation 11: 



^ /./.(*. id] =■ 



V ; J 



[n 7 ,/„(*,n)]=- 



[dcp j 

and likewise for the (easier) case of f c . 

Somewhat more novel is the fact that: 



(15) 



(16) 



ai w((p,x) = [a,,exp(£(<p. +in .)d) -|(...))]|0) + exp(...)a ! |0 
= {(p i +iK i )w((p,K) + 



(17) 
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Thus: 

Trif&t) = \Tr(p ai + pal) = \ < w" a t w> +\< (a^)" w > 
= \ < <Pi + in f >+\<(p i - in i >=< <p. > . 

and similarly for IT and 7ij. More generally, we may deduce that: 
Tr(pf n (®,n))=<f((P,x)> (19) 
for any analytic function f. Intuitively, the effect of the normal product is to keep all the 
aj + terms to the left, where they yield factors of ((pi-i7tO, and to keep the ai terms to the 
right, where they yield factors of 

Finally, consider the classical dynamical equations for the system defined by 
equation 1: 

*,=-¥- (21) 

dqtj 

Taking the expectation value of equation 20, and exploiting equations 15 and 19, we may 
deduce: 

< <j>. >= -2/Tr(p[<I> j,H n (®,Il)] (22) 
Likewise: 

<n j >=-2/Tr(p[n.,H„(^,n)] (23) 
It is tempting to simply interpret equation 22 to mean: 

"4>."=-2z [$.,//] (24) 
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and declare that we are all done (aside from scalar factors which, like h=c=l units, do not 
change the physics.) Note the implicit assumption - as in real- world QFT 1 ' 4 - that the 
Hamiltonian operator H actually used by QFT is the normal-product Hamiltonian. 

To be more rigorous, one must follow the procedures given by Weinberg 1 . The 
operators 4>j and IT above correspond to the q n (x,0) and p n (x,0) operators in his 
development. The true Heisenberg operators, Q n and P n ,are simply defined by his 
equations 7.1.28 and 7.1.29. Equations 22-24 basically tell us that Q and P - which obey 
the Heisenberg-Liouville equations automatically - also fit the dynamics at any time 
"0", and thus that Tr(p(t=0)Qj(t)) yields <cpj(t)>, using the exact same similarity 
transform argument that Weinberg invokes. 

Notice that we only require that Qj(t), so defined, has the property that 
Tr(p(0)Qj(t))=<cpi(t)> for any physically realizable matrix p(0). There is no reason why 
Qj(t) should have to be unique in having that property. In fact, it is possible to derive 
alternate dynamics for an alternate Qj(t) matrix which possesses this same property. This 
is related to the concept of "null equivalence," which is briefly introduced in a previous 
paper 2 , but not central to the goals of this letter. 

IV. IMPLICATIONS AND POSSIBLE EXTENSIONS 

These results indicate that there is an exact equivalence between the dynamics of a 
bo sonic QFT, in the 0+1 -D case, and the statistical dynamics of the corresponding 
"classical theory," with some minor adjustment of scalar parameters. The extension to 
3+1-D should not be difficult; see Appendix B. 
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Does this imply that bosonic QFTs are totally equivalent to classical field theory? 
Bell's Theorem experiments tell us that a complete equivalence is impossible 2 . However, 
QFT is made up of two major components, two bodies of assumptions which are 
combined together to generate predictions - quantum dynamics and quantum 
measurement. Classical field theory was not so explicit in breaking out assumptions 
about dynamics versus assumptions about measurement, but it, too, made strong apriori 
assumptions about the role of causality and statistics in measurement which were not 
derived at all from the assumed Lagrange-Euler equations. Therefore, it makes sense to 
conclude that the differences between QFT and CFT are mainly due to differences in 
measurement formalisms and assumptions about measurement, rather than differences in 
dynamics. Many authors have argued that quantum measurement could in fact be derived 
somehow from quantum dynamics; if true, this might imply that Lagrangian PDE models, 
governed by Heisenberg dynamics, would actually obey quantum measurement rules as 
well. 

The classical-QFT equivalence in dynamics probably extends well beyond 
bosonic field theory as such. LeGuillou et al (hep-th/9602017) claim to provide a general 
procedure for mapping fermionic QFTs into equivalent bosonic QFTs, based on an 
extension of a general procedure published by Witten in 1984, which they cite. This 
suggests that the dynamics of the standard model itself may be equivalent to a bosonic 
QFT, which in turn is equivalent to a classical model. 

This then suggests a possible strategy for developing a finite (not only 
renormalizable) variation of the standard model, by exploiting the classical-QFT 
equivalence in dynamics. 
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Many of the renormalizations required in QFT are based on effects like the 
infinite self-repulsion of a point charge. The same problems exist in the corresponding 
classical theories. But in classical field theories such effects can be eliminated by 
representing particles as extended bodies - more precisely, as true topological solitons . 
In the electroweak part of the standard model, particles mainly derive their mass from the 
Higgs field, a bosonic field whose detailed properties are still almost unknown 6 ; it would 
be reasonable, then, to modify the Higgs model by adding a Skyrme term, which would 
"explain why massive particles exist" (something closely related to explaining the mass!), 
and would ensure the good behavior of the corresponding classical PDE in tasks such as 
the prediction of scattering amplitudes. The PDE-QFT equivalence could then be 
exploited to prove the good behavior of the QFT. (A few auxiliary terms, analogous to 
Fadeev-Popov or Umezawa terms, may be needed to strictly enforce the topological 
constraints assumed in such models.) 

This is only a possible direction for future research; it is certainly not a proof! 
However, please recall that when Weinberg and Salam first proposed the electroweak 
model in the 1960s, they did not have a proof of renormalizability - only a plausible basis 
for hope. The standard model is renormalizable, despite the pessimistic earlier predictions 
based on power-counting rules of thumb, because special properties of the model 
(symmetry) allow good behavior even when power-counting rules are not met. It is 
reasonable to hope for a similar sort of special situation here as well. 

The discussion so far in this section only addresses EWT. But there is no reason 
to believe it could not be extended to strong nuclear interactions, the domain where 
Skyrme models have proven most useful in practical applications. (Witten and others 
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presented many intriguing suggestions for how to develop more empirically-grounded 
models of strong interactions, using topological soliton models, in Chodos et al. 7 ) It may 
even be possible to modify such a finite standard model, by applying the same methods 
John Wheeler used to convert Maxwell's Laws into an "already unified field theory." 
Thus there may even be some hope, in the long-term, of achieving finiteness and 
unification in a field theory which does not ask us to assume the existence of additional, 
speculative dimensions of space-time. 
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APPENDIX A: RELATION TO METHODS IN QUANT-PH 0008036 

The classical density matrix p defined above is related to some of the representations 
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given in quant-ph 0008036, but is quite distinct. The relationship is analogous to 
the relationship between Fourier analysis and Laplace analysis. Fourier analysis, 
Laplace analysis and classical perturbation theory are not alternative theories of the 
universe; rather, they are alternative methods for analyzing the behavior of ODE and 
PDE; they are all sometimes valid and sometimes not, even as applied to well-posed 
ODE and PDE, depending on the properties of the specific ODE and PDE. There is no 
apriori reason to expect that every well-posed model of nature must fit any one particular 
choice of method, such as perturbation analysis; indeed, the great success of 
nonperturbative methods in recent years demonstrates the value of being open to a variety 
of methods. 

The paper quant-ph 0008036 provides a number of alternative representations or methods 
which may appear quite complicated at first. This appendix will summarize the 
two representations used in that paper (v and z) which are closest in spirit to p. 

The matrix p has been defined here as an object over the usual Fock-Hilbert space 
of bosonic QFT. That space is essentially just the vector space spanned by the vacuum 
state I0> and by all states which result when we multiply that state by the usual 
creation and annihilation operators, ai + , ... a n + ,ai,...a„ in any sequence. However, instead of 
representing the statistical moments as a matrix over the conventional Fock space, it is 
possible to represent them as a vector over an expanded Fock space. We can define new 
creation and annihilation operators, bi + ,...,b n + and bi,...,b n , which commute with the usual 
system of ai + ,...,a n . The expanded Fock space is the space spanned by the states we get 
when we multiply any of the "a" or "b" vectors by the vacuum state, in any order. 



Let us consider how this works out when we analyze the statistics of a second- 
order ODE system defined by: 

n 

L(p,0) = |0 2 -il»j 2 -f(<P) . (25) 

7=1 

where f is an analytic function. In CFT, this directly yields the Hamiltonian and 
Lagrange-Euler equations: 

H(^,0) = i0 2 +i|», 2 +/(?) (26) 

7=1 

<p ; . =-w>,- -3, ■/(£), 7 = 1 /i (27) 

This is actually far more general than it seems, insofar as we can use the components 
of $ to represent field values at points in some kind of periodic lattice. Note that we have 
the freedom to pick the Wj parameters to be anything we like, and modify f accordingly, 
to arrive at equivalent representations of the same system. Likewise, we can use this to 
represent either locations or frequency modes in the periodic crystal. 

The statistical moments of this system may be represented by use of the following 
vector over the expanded Fock space, for any pure state: 



v(<p,<p) = exp 



f n \ 

+<Pjbj) 

V ^ J 



1 > (28) 



A statistical ensemble or probability distribution Pr((|>, 3 t $) may be represented by: 

v = $v((p,(p)Pv((p,(p)d n (pd n (p (29) 

In order to simply or "reify" the analysis, we may define another vector z over the 
expanded Fock space: 

z((p,(p) = Xv((p,(p) , (30) 
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where X is the "continuous reification operator" defined by: 

f _ n \ 



X = exp 



— V (a ,a , + a] a] + b b , + b^b*) 

g y j j j j j j 11' 



(31) 



This concept of reification was described at length in quant-ph 0008036, and some of the 
prior papers it cites. Equation 31 is somewhat more convenient than the specific forms of 
reification chosen there, for several reasons. This operator has the following really critical 
properties, which can be verified easily by brute-force algebra, expanding equation 31 as 
a power series and using well-known properties of commutators: 

XaiX" 1 = ( ai + a0/sqrt(2) (32) 

X ai + X" 1 = ( ai + - ai )/sqrt(2) (33) 
The same equations also hold if "a" is replaced everywhere by "b". Again, see quant-ph 
0008036 for generalization to the PDE case. 

Applying well-known methods of commutator algebra to equation 28, we 
easily deduce: 

aj v ($, a t $) = wj (pj v($, 3 t $) (34) 

bj v ($, 3 t $) = (3 t (pj) v($, a t $) (35) 
and, in general: 

/({w,<pJ,<p)vX<p,<p) = f{a,b)v(q,<p) (36) 

In this case, the definition of f(a, b) is unambiguous; it does not require us to specify 
whether we are using the canonical or normal or other form, because these annihilation 
operators all commute with each other. 

We can derive the dynamics of v simply by differentiating equation 28 with 
respect to time, which yields: 
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v((p,(p) = d t 



£(ir,(p,a; +(p,b[ exp Y,( w j<P j a ) + <P j h ) 



I0> 



= L +<p j b])v((p,(p) 

j 

Exploiting equation 27 and defining 

gj($) = -djf($), 
equation 39 reduces to: 

v(y,<p) = Y,{w j (p j a) +b + j (-w)(p j +g 7 (9)))v(9,£) 



(39) 



(40) 



(41) 



If we substitute in from equations 34 and 35, and then exploit linearity to sum over all 
states in a statistical ensemble, we may deduce for any physically realizable vector v that: 
3 t v = Go, v v + Gi, v v (42) 



where: 



G o,v =Y, w M + i b i - h : (t : » 



(43) 



and: 



G /.v =Y. b Ui 



(I w 
a 



(44) 



where I apologize for the awkward notation in equation 44, where the angle brackets 
are used to denote the n-component vector made up of the operators (aj/wj). 

To derive the dynamics of z, for any pure state, we simply insert the definition in 
equation 30 and exploit the relations in equations 32 to deduce: 
3 t z = Go, z z + Gi, z z , (45) 



where: 



Gq,z - Go, v 



(46) 



14 



and 

(47) 



where I now avoid awkward notation by defining the operator vector 4> as the vector 
made up of the n components: 

a, + a* 

<$>j = 1 J. (48) 

Note that the operator G z is antiHermitian. Thus the vector z obeys a Schrodinger 
equation with all the standard properties assumed in most theoretical discussions of 
quantum systems. If quantum measurement could be derived as a consequence of 
quantum dynamics, in principle, then any classical system in this class should actually be 
governed by quantum measurement as well; however, in recent discussions with Shih, 
Dolmatova, Dowling and Strekalov, we have found some evidence suggesting that the 
real story with quantum measurement is even more complicated than the story given in 
section 2 of quant-ph 0008036 - more complicated for any dynamical system, classical or 
Heisenberg or Schrodinger. 

In the v and ^representations, the energy operators are not the same as the "gain" 
operators G v and G z . Let E represent the classical scalar energy of a system as defined by 
equation 26. 

There are two types of energy operator for statistical representations like v and z. 
I will call these X-type energy operators and Q-type energy operators. 
H is a X-type energy operator for v if and only if: 

H v = E v (49) 
for all pure states v. H is a Q-type energy operator for v if and only if: 
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v H H v = E Ivl 2 



(50) 



for all pure states v. The definitions for z are the same, with "v replaced by "z". 
It is easy to see that all X-type energy operators are also Q-type operators, but that 
the converse is not in general true. 

From equations 26 and 36, it is easy to deduce that equation 49 is satisfied by the 
operator: 



where again I apologize for the awkward notation. As before, equation 31 effectively 
induces a similarity transformation on equation 51, such that we may deduce that the 
following operator is a X-type energy operator for z: 



Note that any statistical ensemble or linear combination of pure states will also be an 
eigenvector of H z , with the same eigenvalue E, if and only if all of the states in the 
ensemble have the same energy E. 

This is clearly getting us very close to the way that QFT calculates spectra!! 
In fact, H z is identical to H Q , except that Hq replaces the leading term on the RHS of 52 
by "7t 2 " term which represents I3 t (|>l 2 in terms of the operators aj and aj + only! If one were 
to apply this formalism to a simple example, like the sine-Gordon system, one would 
observe that the kink soliton, in the rest state, has a G z eigenvalue of zero (is not 
oscillating!) even though it has an H z eigenvalue corresponding exactly to its classical 
energy. 




(51) 



(52) 
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To conclude this appendix, I should add some additional comments on the 
concept of null-equivalence which was introduced briefly but not exhaustively explored 
in quant-ph 0008036. That paper introduced some alternative density matrices, such as 
p z , which are quite different from the p defined above. The complete analysis requires 
using a concept of "striated matrix" - essentially, a matrix where elements like pij match 
elements like Pi+j ;0 as required for any physically realizable matrix. Any proposed density 
matrix is null-equivalent (as defined there) to one and only one striated matrix, which 
may or may not be physically realizable. A striated matrix is physically realizable as a 
pure state if and only if it is rank one. When an operator is null-equivalent to a valid A,- 
type energy operator, it is itself a valid Q-type energy operator. Again, these comments 
are intended to fill in the analysis of v and z representations given in quant-ph 0008036, 
and should not be taken out of context. 

APPENDIX B. CLASSIC AL-QFT EQUIVALENCE FOR BOSONIC QFT IN d+1 
DIMENSIONS 

Quant-ph 0008036 provided excruciating detail for extensions of the v and z 
and other representations to the case of d+1 dimensional systems. This appendix 
will discuss a few key details of how to extend the derivation of Heisenberg dynamics in 
section III of this paper to the d+1 dimensional case. First let us try to analyze the 
statistical dynamics of the following class of CFT, expressed in Hamiltonian form: 

L = \Y^{(j> j K j -TijVj-Wpj I 2 -m 2 <p 2 )-/(p,;r,Vp), (53) 

7=1 

where $ e R n and f is an analytic function. In theoretical physics, we usually restrict 
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our attention to the special case where the mathematical vector $ is actually composed of 
relativistic scalars, vectors, iso vectors, tensors, etc.; however it is easier here to address 
the general case. Whenever the original Lagrangian is relativistically invariant, then the 
equivalent QFT will also yield predictions consistent with relativity. 

Note that equation 2 has some redundancy in it. Any dynamical system with a 
Lagrangian of this form can be expressed in many equivalent ways; for example, we can 
change the "mass" terms mj, and balance those changes by changes in f, to arrive at 
a "different" Lagrangian which is actually the same functional of the fields. (In principle, 
we could even make the mass terms functions of Igl, after a Fourier transform; however, 
I have no need of that here.) There is no one "right" choice between these alternatives; 
all are valid, so long as they are consistent with later stages of analysis. 

CFT then gives us the classical Hamiltonian density and Lagrange-Euler 
equations: 

H =i£(j V( Pj I 2 +m 2 J <p))+f(<p,n,V<p) (54) 

7=1 

*>-w. <55) 



Ttj =A(p j -m 2 j (p j - 



f 



Sf 



(56) 



Let S(t) denote the state of this dynamical system at any time t; S consists 
of all of the values of $(x,t) and 7t(x,t) across all points in space x at time t. Then 
the classical density matrix may be defined by the functional integral: 

C v(S)v"(S) 

p= — ~ 7 — Pr(5)J S, (57) 
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where: 



v(S) = exp 



o , 



(58) 



where d is the number of spatial dimensions (i.e. xeR) and: 



9 ] (p) = ^^)\e- il - y -(p j (y)d d y_ 

t j(p^ = ~r== \ e '- ln j(y) dd y 



w j (p) = Jm j + 



(59) 
(60) 

(61) 
(62) 



It is straightforward but tedious to then prove (for d=3) that: 

Tr(p^,(x))=<<p,(x)> (63) 

Tr(pn 7 (x))=<^,(x)> (64) 

where the field operators <£j and Ilj are defined exactly as in Weinberg, equations 
7.2.29 and 7.2.30, which in my notation may be written: 

n J (x) = n + J (x)+n-(x) 

where: 

^- j (x) = {^ + j (x)f (67) 

n-(x) = (n + J (x)) H (68) 



(65) 



(66) 
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(69) 



U + j(x) = -ic J" {jwj(p))e '--a j (p)d d p 



(70) 



In actuality, to perform these calculations, I first treated the parameter c in equation c 
as an unknown parameter, and solved for ^j + (x)v(S)=((Pj(x)+ imaginary)v(S), 
and likewise for Tl. (I found it convenient as well to replace "p" by "q" in equations 58- 
61, before proceeding. Those definitions could be expressed a bit more elegantly, but the 
form given here was convenient for derivation.) 

With these definitions, one would expect the logic behind equation 19 to go 
through as before, yielding the exact same equation for the PDE case! The standard 
commutator relations from equations 10-16 would also go through (adjusted to the PDE 
case) - but may simply be extracted from standard texts like Weinberg, since these are 
the standard field operators now. An important example is: 



where once again the "n" subscript refers to the normal product version of the function. 
The derivation should then go through, without the previous caveats about needing to 
adjust scalars. 



[<*> j (x ), { /„ my),U(y))d d y] = id d (x - y) 1-$- my), U(y) 




n 



(71) 
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